Response of nucleons to external probes in hedgehog models: 
II. General formalism 
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Linear response theory for SU(2) hedgehog soliton models is developed in analogy to a standard 
method in many-body physics. In this framework, we discuss response of baryons to external probes, 
and develop expressions for polarizabilities. We discuss isospin effects (neutron-proton splitting) in 
polarizabilities. Methods for cases with zero modes are presented, including numerical techniques. 
Our approach is based on the l/JV c -expansion scheme. We work in a model with quark and meson 
degrees of freedom, but the basic method is valid in any hedgehog model, such as the Skyrmion 
or the Nambu-Jona-Lasinio model in the solitonic treatment. The equations of motion for coupled 
RPA quark-meson fluctuations are classified according the hedgehog symmetries, and are written 
down explicitly in the grand-spin basis. 
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> ! I. INTRODUCTION 

m 

■ In recent years various hedgehog models (Chiral 
Quark Meson (CQM) models @|f §@ lli , Skyrme mod- 
els @,|,|,|o|, hybrid _bag models chiral mod- 

2 yi3|fl4[ ] or the Nambu-Jona- 
151 in the solitonic treatment 



C^l ' els with confinement L 
Lasinio (NJL) model 

@0JT|Jl|,|2^,|l|,|2|) were extensively applied to de- 
Oh! scribe the physics of low-energy baryons. Semiclassical 
(JL) ■ methods for treatment of these models, such as various 
qj[ projection methods H^|,||,|§, or RPA method §7]] 
were developed. Masses, various charges, tt — N phase 
shifts |2q] , were calculated, with quite reasonable agree- 
ment with experiment, depending on the specific model, 
number of included fields, etc. In this article we develop 
the linear response formalism for hedgehog models. We 
work in the framework of a CQM model, since it has 
both quark and meson degrees of freedom, and in this re- 
spect has the essential features of both the purely mesonic 
Skyrme model, and the NJL model, which involves quark 
degrees of freedom only. Our methods and final expres- 
sions can be modified straightforwardly to be applicable 
in these models. 

Hedgehog models can be used to describe response of 
nucleons to external probes, and to calculate correspond- 
ing polarizabilities. A natural approach is the linear re- 
sponse method of many-body physics |29]. The under- 
lying picture is as follows: A current interacts with the 
nucleon, creates an intermediate state which is an RPA 
phonon excitation on top of the soliton. This state in- 
teracts with another current and de-excites back into a 
nucleon state. The RPA phonon states are constructed 
from one-particle-one-hole excitations of the quarks, as 
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well as from quantum meson excitations. Quark and me- 
son fluctuations are coupled, and the resulting equations 
of motion for the fluctuations are solved. An example of 
physically important two-current observables which can 
be calculated in this way are the electromagnetic polar- 
izabilities of the nucleon ||(| . This topic is extensively 
studied in the preceding paper (3^] , henceforth referred as 
(I). The present article is devoted to development of the 
necessary formalism, and contains many technical but 
necessary details of linear response in hedgehog models, 
while (I) concentrates on physical aspects. 

This article is organized as follows: In Sec. || we very 
briefl y rev iew a CQM model |3^], its soliton sol ution s 
(Sec. II A), as well as hedgehog symmetries (Sec. II B). 
One of the discrete symmetries, the grand-reversal sym- 
metry |35) , |23| ], will be particul arly useful in classifying 
various perturbations. Section III is the core of the pa- 
per, and describes the equations-of-motion approach H] 
to linear response in hedgehog models. We start from dc- 



riving small-fluctuation equations of motion (Sec. [II A) 
for coupled quark-meson systems driven by an external 
perturbation. These equations are classified according to 
hedgehog symmetries. In the static limit the grand re- 
versal symmetry, 1Z, decouples the equations into odd-7?. 
equations, involving quarks only, and even-7?. equations, 
involving both quarks and mesons. We d iscuss in de- 
tail the problem of zero modes (Sec. [II B). These zero 
modes arise from braking of the symmetries of the la- 
grangian by the soliton solution. In our applications, we 
will have to deal with rotational zero modes (cranking) 
and translational zero modes (isoscalar electric pertur- 
bation in (I)). We describe a numerical method to deal 
with the excitation of zero modes whose amplitude di- 
verges as the frequency of the perturbation goes to zero 
(Sec. IIIC ). We discuss stability of solitons (Sec. IIID| ). 
In Sec. [HE we present cranking in the linear response 
form alism. Q uantiz ation via cranking is reviewed in Sec. 
Ill F . In Sec. IIIG we describe the calculation of polar- 
izabilities in states of good spin and isospin, and obtain 
our basic formulas. Section IV illustrates the method 
by presenting the standard calculation of the N-A mass 
splitting, as well as the evaluation of the neutron-proton 
hadronic mass difference. The issues of A c -counting are 
discussed in Sec. |y| . We show how to apply the linear re- 
sponse in a way consistent with l/iV c -expansion scheme. 
Finally, Sec. |v| contains remarks relevant to other mod- 
els (Skyrme model, NJL). 

Appendices contain some details of the grand-spin al- 
gebra, derivation of the explicit forms of the equations of 
motion for fluctuations in of the grand-spin basis (App. 
A) , and a glossary of useful formulas with collective ma- 
trix elements (App. B). We also give a simple proof of 
equality of the soliton mass and the inertial mass param- 
eter (App. C), and discuss the issue of Pauli blocking of 
the Dirac sea in chiral quark models (Sec. D). 
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II. HEDGEHOG MODELS 



In this paper most of the derivations will be done in 
the framework of the chiral quark-meson model (CQM) 
of Ref. [0. For the details, description of solutions, and 
the resulting phenomenology obtained with the cranking 
projection method, the reader is referred to Ref. [Ejjl . 
The reason of choosing this particular model over other 
models, e.g. the Skyrmion or the NJL model, is that it 
contains both quark and meson degrees of freedom, and 
formally has all essential features of a generic hedgehog 
model with two flavors. At the same time, it is free of 
the non-linear complications of the Skyrme model, or the 
Dirac-sea complications of the NJL model. 



A. Soliton solutions 



The lagrangian of the model is the Gell-Mann-Levy 
lagrangian pq] , with ip denoting the quark operator, and 
a and 7r denoting the meson fields: 

C = [i<j) + g (a + ij 5 t ■ 7r)] if) 

+ i(^<r) 2 + |(^7T) 2 -U(a,ir). (2.1) 

The Mexican Hat potential, 




leads to the spontaneous breaking of the chiral symmetry 
in the usual way ||3JlJ]. Our convention for the pion 
decay constant is = 93MeV. At the (time-dependent) 
mean-field level, only valence quarks, denoted by q, are 
retained in the expansion of the quark fields, and the 
meson fields are treated as classical, c-number fields p^] 
(see also App. D). The time-dependent equations of 
motion have the form 



(/#] - id t )q = 0, (2.3) 
-a<p=^M-gN c qMq, (2.4) 

where cf> — (er, 7r) denotes the meson fields, M — (/3, «75t) 
describes the quark-meson coupling, and the Dirac hamil- 



tonian is h[4>] = —ia ■ V — gMtfi. Equations (2.4) have a 
stationary solution of the form 

a(r,t) = a h (r), n(r,t) =rn h (r), q(t) = q h {r)e~ le \ 

*=(i CT G fi?(r)) (|ui) - |rfT))/ ^' (2 ' 5) 

where e is the quark eigenvalue. For discussion of this 
solution, plots of the radial functions ah, ^h, Ght and Fh, 
and other details, the reader is referred to Refs. p],B3||. 
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B. Hedgehog symmetries 



The solution (2.5) has the hedgehog form, which breaks 
the spin, J, and isospin, J, symmetries of the lagrangian 
( |2.lD , leaving as a good symmetry the grand spin, K — 
I + J. There are also two discrete symmetries which 
are very useful in classifying solutions and perturbations. 
One is parity, T 3 , the other is the "grand-reversal" sym- 
metry, TZ, discussed in Refs. [|HJ2^]. Formally, TZ is de- 
fined as the time-reversal, followed by an isorotation by 
angle ir about the 2-axis in isospin. Explicitly, it trans- 
forms the quark spinors and mean meson fields as follows 

q(r,t) -> a z T 2 q*(r,-t), 

a{r,t) -> <r*(r,-t), n{r,t) w*(r,-t). (2.6) 

We denote the action of TZ on an object by the super- 
script n . The soliton solution has K vn = ++ . 



III. LINEAR RESPONSE IN HEDGEHOG 
MODELS 

In this section the basic formalism of linear response in 
hedgehog models is developed. We use the equation-of- 
motion approach |34f| , which is based on solving equations 
of motion for small oscillation on top of the ground state 
solutions. This method is equivalent to the particle-hole 
formalism p9| , in which one introduces a quantum RPA 
state, quasi-boson RPA phonon operators, etc. Meth- 
ods such as cranking, projection, or quantization of zero 
modes, can be described in this framework, and have def- 
inite quantum-mechanical interpretation. For simplicity 
of notation, we present our formalism in the equations- 
of-motion method. 



A. Equations of motion for small fluctuations 

Let us consider a small oscillation problem in our sys- 
tem. We introduce shifts in the valence quark spinor and 
in the meson fields, 

5q{r,t) = (X(r)e- iuJt + Y n {r)e iu>t ) e~ iet , 
5<f> a (r, t) = Z a {r)e~ iut + Zf{r)e^\ (3.1) 

where X and Y describe the shift in the valence quark 
spinor, and S(f>o and Scf> are the shifts in the a and 7r 



fields, respectively. Note, that in Eqs. (3.1) the TZ trans- 
formation has taken the place of the usual |^9| complex 
conjugation. This is because in hedgehog systems the 
grand-reversal replaces the u sual time-reversal symme- 



try. According to definition (2^), the meson shifts 5<j) 
are even under grand-reversal, but the quark shifts have 
in general both even and odd components. We linearize 



equations (2.4) about the solitonic solution (2.5), and 



obtain the quark- meson RPA equations. When external 
perturbations are present, these equations are in general 
driven by a quark source, j q , and a meson source, 



(3.2) 



Again, the meson source is even under TZ, whereas the 
quark source has in general even and odd- TZ components. 
Using the fact that h[4>h], M and qt are even under TZ (in 
fact they are K vn = ++ objects), we obtain a general 
form of the linear response equations for our hedgehog 
system: 

(h[<M -e)X-gY J M a q h Z a - uX = jx , 

a 

(h[<f> h \ -e)Y-gJ2 M a q h Z a + ujY = j Y , 

a 

s 2 u 



gN c (q f h MaX + Y^M a q h ) - uj 2 Z a = j z . (3.3) 



Introducing auxiliary meson mom entum variables P a — 
—iujZ ai we observe that Eqs. (3.3) can be written in the 
symplectic form E9] 



HZ - cjA£ = j, 



(3.4) 



where H is the RPA hamiltonian, and A is the symplectic 
RPA metric, satisfying A 2 = 1. In th e gr and-spin basis 
(App. A), 7i is real. Our problem (3.3) can then be 
written as 



n 



( N c {h-e 

V 



-gN c Mq h 0\ 

N c (h-e) -gN c Mq h 

gNj h M -gN c q\M -V 2 + U" 

1 / 



(3.5) 



A : 



10 
0-100 
i 
-i 




(3.6) 



Note the appearance of an odd- TZ momentum component 
in the source, jp, which arises in some cases (cranking). 

We are interested in the limit of vanishing frequency 
of the external perturbation, lu — > 0. If zero modes are 
excited by an even- TZ perturbation (Sec. IIIB), then the 
full equations Eq. (3.3) have to be solved. Otherwise, 
one can set uj — and deal with simplified cases. At 
this point the grand-reversal classification becomes very 
useful. Acting with TZ on Eq. (3.4) effectively replaces 
X 



Jp- 



Y, jx <-> jy, Z <-> Z, j z <-> jz, P- P, and 
-> jp . Let us introduce odd and even grand-reversal 
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comb inat ions: 5q = X ± Y, jr 1 = jx ± Jy, and rewrite 



Eq. (^4) by adding and subtracting the first two equa- 



tions. We get for the case of an odd-TZ perturbation 

(h-e)Sq~ =j~, 

P = jp, (3.7) 

and for the case of an even- 1Z perturbation 

(h-e)5q+ -2gMq h Z = j+, 
(-V 2 + U")Z - N c gq) i M8q+ = j z (3.8) 



The odd- 1Z equations ( |3.7[ ) involve a quark field equa tion , 
and a trivial equation for P. The even- 1Z equations ( |3.§| ) 
involve coupled quark and meson fluctuations. Equations 



(3.4), or (3.7,3.8) are further decomposed by grand-spin, 



K, and parity, V (App. A). 

In models with vector mesons, such as |32 , the odd- 1Z 
equations may also involve mesonic shifts. For example, 
the space components of the w meson and the time com- 
ponenent of the p meson enter into the cranking equa- 
tions of motion [|36|]. 



B. Zero modes 



First consider the undriven problem (3.4), with j = 0, 
which determines the RPA spectrum and eigenmodes. A 
complication arises whenever a continuous symmetry of 
the lagrangian is broken by the solitonic solution, e.g. 
translation, or rotational symmetry. For each broken 
symmetry the small fluctuation equations have a pair of 
zero-modes p9fl : £o, the symmetry mode, obtained by 
acting with a symmetry generator on the solitonic solu- 
tion, and a conjugate zero mode, £i. They satisfy the 
equations 

HZo = 0, (3.9) 
H£i = -iA&. (3.10) 

The remaining "physical" modes, £j, satisfy the equations 

W& = WiA&. (3.11) 

One can easily show that the symplectic norms satisfy 
conditions 

^ = ^ = 0, i = 2,3,... (3.12) 

where M. is the appropriate inertia parameter (mass, mo- 
ment of inertia) parameter, and M% are the symplectic 
norms of the physical modes. The factors of j are con- 
ventional, and factors of i are inserted for convenience. 



Expanding the solution of Eq. (3.4) in RPA eigenmodes, 



G 



£= E (3-13) 

/i=0,l,2,... 

introducing "charges": Q = 4i^j, = 4£tj, // = 
1,2,..., and using Eqs. ( [3.12] ), we find that 

ujc\M + Qo = 0, iujcqM + Qi — c\M = 0, 
d{ui - = Qi. (3.14) 

We consider two cases which arise in practical applica- 
tions: 1) Qi = 0, and 2) Q x ^ 0, Q = 0. 



1. Case 
Using Eq. (|3.14|) we find 



«<3o Qo Qi fo -i r\ 

.Mw 2 Ali(uJi — oj) 

The second-order energy shift, k, corresponding to a 
given perturbation (a "polarizability" is equal to 2k) is 
given by the usual perturbation theory result 



ti- 
i Qo 



2 Mw 2 ' " 2 ^M{^-^Y 

In the limit to —> 0, the coefficients Co, ci and the zero- 
mode part of K diverge, as long as Qo ^ 0. This has a 
physical interpretation: for instance in the case of trans- 
lation the center of mass of the system moves, and the 
amplitude of this motion, Co, as well as "velocity", ci, 
diverge. In (I) we show how this feature of the linear 
response formalism leads to the Thompson limit of the 
Compton scattering amplitude. 



zero , phys. 
= K + K , 



2. Case Q 1 / 0, Qo = 



In this case we ca n take the limit u> — > on the outset, 
and from Eq. ( 3.14 ) we get 



ci = 



Qi 



Q, 



Mi 



(3.17) 



The amplitude of the symmetry mode, Co, remains un- 
dermined. The second-order energy shift is: 



l Qi , l \ "* 

2 M 2 ^MiWi 



Ql 



(3.18) 
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C. Numerical methods in presence of diverging zero 
modes 



Numerically , the ex citation of amplitude-growing zero 



modes (Sec. [II Bl) creates special difficulties in ex- 
tracting the "physical" parts of observables, e.g. elec- 
tromagnetic polarizabilities. The pro blem can be reme- 



died as follows: We solve Eqs. (3.4) for a small value 
of lu. Next, we project out the zero-mode part from £, 
obtaining ^ h y s - = £ — co£oj and calculate physical parts 
of observables. The procedure is repeated with decreas- 
ing u, until the results no longer change. In practice, 
a very high accuracy of the soliton solution, as well as 
the fluctuation solutions, is required for this procedure 
to be feasable. A better method is to project the part of 
the source, j, which couples to the zero mode, and solve 
equations 

nS Phys. _ ^^Phys. = jphys. ^ (3 1Q ) 

where jP h v- — j — (Q /M)A£i, an d gi is obtained by 
solving Eq. (3.10) first. Equations ( |3.19 ) do not excite 



the zero mode, and directly lead to the physical part 
of the solution. The advantage of the method with the 
projected source over the direct method described previ- 
ously follows from the fact that in numerical solutions to 



Eqs. (3.1E) the admixtures of the zero mode arise only 
from numerical noise. Their amplitude is small, such that 
we can easily control numerical precision in the physi- 
cal mode. Because of these admixtures, a small nonzero 



value of u should be kept as a regulator in Eqs. ( p.l9[ ), 
and the zero-mode contamination has to be projected out 
after the numerical solution is found. 



D. Stability of solitons 



Since in our problem TL and A are hermitian, one finds 
that = ujfii is a hermitian eigenvalue problem. 

Therefore in our case u>f are real, and u>i can either be 
purely real, or purely imaginary. The modes appear in 
conjugated pairs with cj,; = — luj. If the spectrum 

contains an imaginary eigenvalue, we have to instability 
(in the Lyapunov sense J37j) of the ground-state (soliton) 
solution pq , p9[ , and of course linear response on top of 
an unstable system makes no sense. In Ref. p7| we have 
shown that the soliton of Ref. jl) is stable with respect to 
breathing modes, i.e. the K p = + excitations. With the 
explicit forms of the equations in App. A, stability could 
be checked numerically for any K p vibrational mode. It 
is generally believed that the hedgehog solitons are in- 
deed stable, although it has not been proved analytically 
or numerically. 
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E. Cranking as linear response 



Cranking 123] may be viewed as linear response. In 
a frame iso-rotating with a small angular velocity A, we 
discover equations of the form ( |3. 10 ) , with u> = and 
j = — iAA£o(A). In this case £o(A) is the symmetry mode 
obtained by acting on the soliton fields with the generator 
of isorotation about the axis A: 



(%/2t- Xq h \ 
i/2r ■ Xqh 
X x n h 




(3.20) 



Next, we have to find the conjugated mode, by solving 
the seco nd of Eqs. ( Off ). We notice, that this is an odd- 



R case fl3.7|). We immediately get P — ^X x Tv h . For the 



quark shift, 5q cr , a differential equation of the form (3.7) 
is sol ved [p3| . The problem is of the type discussed in 
Sec. [II B 2 , where M. is the moment of inertia, 0, and 
the "charges" are: Q 1 = A6, = for [i ^ 1. The 
second-order energy shift is: K = ^A 2 0. Explicitly, one 
finds 



e = e m + e„ 

&m = J d 3 x(X x Tv h ) 2 = (8tt/3) J drr 2 TT 2 h , 
Q q = 2 / d 3 xSqt r X ■ rq h (3.21) 



F. Quantization 

The simplest approach to quantization via cranking, is 
to recognize that in the frame isorotating with velocity A, 
in wh ich we solve the cranking equations of motion(Sec. 



HIE ), we still have the freedom of (iso)rotating the soli- 
ton by an arbitrary (time-independent) angle. This is 
an example of the freedom of choice in the cq coefficient 
in Sec. |III B % which in this case corresponds to three 
Euler angles, or, in the commonly used Cayley-Klein no- 



tation pi, to the matrix B = bo + b ■ t |23 . In our 
mean-field approach, the corresponding fields carry these 
(time-independent) B matrices, and in the rotating frame 
they assume the form: 



(T, 



BtzB\ q^Bq. (3.22) 



Matrix B plays the role of coordinate variables conju- 
gated to A, which upon quantization becomes a differ- 
ential operator P,p3[. The quantization is straightfor- 
wardly implemented in two steps: 1) one identifies the 
collective spin and isospin operators, as done in Ref. p3[ . 
Then 

A9 = J, I a = c ab J b , (3.23) 
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where J and / are the spin and isospin operators, sat- 
isfying appropriate commutation relations, and c a b, de- 
fined in App. B, has the meaning of the transformation 
matrix from the body-fixed to the lab frame [^9). 2) 
Corresponding collective wave functions are introduced. 
Expectation values of operators are calculated by first 
identifying in the semiclassical expression for an opera- 
tor its collective part (dependent on A, c a b, etc.), and an 
intrinsic part (dependent on the meson and quark fields 
a, 7r, q). Then, the matrix element factorizes into a col- 
lective matrix element in the wave functions of App. B 
(this is an integral over the collective coordinates, viz. 
Euler angles, or (bo, £>)), and an intrinsic matrix element, 
which is a space inte gral over the quark and meson fields. 
For details, see Ref. P3|| . 



G. External perturbations 



The quark and meson field profiles in Eq. ( 3.22 ) are 
in general not equal to th e hedgehog profiles. We have 
demonstrated in Sec. [II E that the quarks develop shifts 
upon cranking. If some other (external) interaction is 
present, then the profiles are additionally shifted. These 
shifts are obtained by solving the linear response equa- 
tions, as described in Sec. III. We introduce a resolvent 



for the TL — toA operato r in Eq. (3.4) (RPA propagator) 



and solve formally Eq. (3.4), obtaining 
t = Qj, Q = (H-loK)-\ 



(3.24) 



In the presence of cranking and some other external per- 
turbation, we have 



£ £,cr ~t~ ^,ext — Q(jcr jext), 



(3.25) 



where subscripts cr and ext refer to cranking, and an ex- 
ternal perturbation, respectively. The second-order en- 
ergy shift corresponding to a perturbation can be written 
as 



- 2£tj = 2 fQj. 



(3.26) 



The difference between this expression, and the generic 
expression (3.16) is that in the present case the source 



jcolljintr^ ThuSj 



carries collective degrees of freedom, j 
the matrix ele ment of k in a baryon state \b) is (see ex- 
ample in Sec. [V A): 



K b = 2(coU\j coU j cM \coll) 

d 3 x d 3 x' f ntr \x)Q(x,x')j mtr (x'), 



2(coU\j coU j coU \coll) d 3 xC ntrf (x)j mtr (x), (3.27) 



where \coll) represents the collective wave function (App. 
B) associated with the baryon state \b). 

It is possible to have isospin-dependent effects in linear 
response of the nucleon. For example, if the external 
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interaction has K v = 1 + (the same quantum numbers as 
in cranking), we pick up cross terms between cra nking , 



and the external perturbation (see example in Sec. tVB): 

K b = 2(coll\j™ u jZ\ l \coll) 

d 3 x d 3 x' C tr Hx)G(x,x')j^l r (x') + h.c. 

= 2(coll\f c ff e ° x l l\coll) J d 3 xCr tr Hx)jT x T(x)+ h.C. 

(3.28) 



Expressions ( [3.27 , 3.28) are just second-order pertur- 
bation results. We may formally continue to higher order 
in perturbation theory, which leads to chains of the form 

K iu ..., in =2(coll\j?f%°^..jtf\coll) 

d 3 X! ... d 3 x n j-" M {x 1 )g(x 1 ,X2)V l2 (x 2 ) 

G{x 2l x 3 ) ... V iri _ l (x n -i)G(x n -i 1 x n )j™ tr (x n ), 

(3.29) 

where Vi k is interaction of k th type. The total energy 
shi ft is t he sum over all possible orderings of 



in ( 3.29 ). Because the groun d state has K v = + , the 



matrix element in Eq. ( p. 29 ) is non-zero only if one can 
compose the K v quantum numbers of Vi 2 , to 



K v = + . In (I) we show an application of Eq. ( 3.29] ) 



with two RPA propagators in the analysis of the neutron- 
proton splitting of electromagnetic polarizabilities. In 
Sec. we discuss in what cases going to a higher order 
in perturbation theory is consistent with iV c -counting, 
which is our basic principle in organizing the perturba- 
tion expansion in hedgehog models. 



IV. SIMPLE EXAMPLES 

In this section we give some simple application of the 
described formalism. A more advanced and physically 
important case of electromagnetic polarizabilities is given 
in (I). 



A. iV-A mass splitting 



As an illustration of application of Eq. ( 3.27 ), consider 
the A^-A mass splitting. In this case n b is the energy shift 
of the bary o n \b) due t o the cranking perturbation. From 



Eqs. (3.23, 3.21, 3.27) we obtain immediately the usual 



expression for the A-A mass splitting: 

M A - M N = i((A|A 2 |A) - (A|A 2 |A))9 = ± (4.1) 
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B. Hadronic p - n mass splitting 



As an example of an isospin-dependent effect, consider 
the neutron-proton mass difference due to the difference 
of the up and down quark masses. The perturbation in 
the lagrangian has the form C m = \(m,d — m u )ipT3ip. It 
has K vn = 1 H , exactly as c ranking, hence a mixed 
perturbation of the form ( |3.28 ) appears. Passing to an 



isorotating frame, we find the source correspo ndin g to the 



quark mass splitting, which arises in Eqs. (3.7): j m — 
\{ m d — w u )A c 7 c • rqfi, where c is defined in App. B. 
Since we have already solved the cranking equation, we 
do not have to solve the new equation with source j m . 
We simply calculate the overlap of j m with the shift in 



the fields due to cranking, Sq cr , according to Eq. (3.28). 
Using the fact that (N\X ■ c|A) = -(N\I 3 \N), we obtain 
the following expression for the hadronic splitting of the 
neutron and proton masses: 



(M n - M p ) hadr - = (m d - m u ) J d 3 x d 3 x'[(n\jlSq cr \n) 

-(n*->p)] 
= m £ _m^ j d s xj inM §q ^r_ (4 2) 

The numerical value, obtained for the solution of Ref. 
H gives (M n - M m ) hadr - = 0.4 x (m d - m u ), which 
for typical values of (md — m u ) gives a number around 
2MeV. The electromagnetic mass difference can also be 
studied in hedgehog models models plifl . 

V. iVc-COUNTING 

The basic organizational principle behind hedgehog 
models is the 1/N C expansion of QCD |p|,p2p3]. In the 



N c — * oo limit, masses of baryons diver ge a s N c , and 
can be calculated using mean-field theory |4^] . It should 
be noted that the assumption of the spin-isospin corre- 
lated wave function, which is essential in hedgehog mod- 
els, does not follow from the large- N c limit alone — it 
is an additional assumption of the hedgehog approach. 
By analogy to nuclear physics, in systems with many nu- 
cleons we may have nuclei with intrinsic deformations, 
but we may also have spherically symmetric nuclei, and 
it is the dynamics which determines whether the wave- 
function is deformed or not. In hedgehog models the 
hedgehog wave function is assumed to be deformed in 
the spin-isospin space, and the nucleon the A masses, 
which are of the order N c , are degenerate in the leading- 
Ac order. 

When cranking is used, these masses split as ~ A^ 1 . 
In fact, cranking becomes an exact projection method 
in the large- A c limit, since it may be viewed as a 
Peierls-Yoccoz projection with (5-function overlaps be- 
tween rotated wave functions Ip29fl. Thus we obtain 



the hedgehog result for the mass splitting, Eq. (4.1) 
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It would not be consistent, however, to conclude that 
the nucleon or A masses individually are given by the 
hedgehog soliton mass plus the cranking piece. There 
are other effects (center-of-mass correction, centrifugal 
stretching, etc.) which enter at the same level as the 
cranking term. Also, the effective lagrangian may be 
supplemented by subleading terms in N c , which we 
did not have to include to obtain the leading piece in 
the hedgehog mass. Therefore, it is useless to write 
down Mj = M h + J{J +l)/{2@) + 0{N~ 1 ), since the 
last term, which we do not calculate, enters at the same 
level as the cranking term. We can only trust the lead- 
ing piece, Mj = Mh + 0(N ( T 1 ), and, in order to main- 
tain consistency with the iV c -counting, the mass for- 
mula should not be "improved" by adding the crank- 
ing term. The mentioned effects of center-of-mass cor- 
rections, centrifugal stretching, etc., are at the leading 
level the same for the nucleon and for the A, there- 
fore for the N — A mass splitting we get the formula 
A/a - M N = 3/(26) + 0{N- 2 ). 

The prescription, which we tried to illustrate above, 
is that with semiclassical methods we can only get the 
leading- N c term for a given observable. The power of N c 
varies, depending on the quantity we are investigating. 
The same is true for the calculation of polarizabilities, 
described in this paper. We can e asily obtain th e N c 
behavior of various terms in Eqs. fl3.27|j3.28|JOs| ), but 



only the leading-iV c piece corresponding to a particular 
polarizability should be retained. As an illustration, con- 
sider the electric polarizability of the nucleon, discussed 
extensively in |5Q| and in (I). The electric field polarizes 
the hedgehog. The electric charge of the quark has an 
isoscalar component, of order N~ , and isovector com- 



ponent, of order 1. We immediately see from Eq. (3.27) 
that the leading part of the electric polarizability of the 
nucleon is obtained from interactions with two isovec- 
tor sources, and the term with two isoscalar sources is 
two powers of N c suppressed. The non-dispersive seagull 
effects also enter at the level of N c (I), hence the nu- 
cleon polarizability goes as N c . Quite analogously to the 
problem of the N — A mass splitting, the neutron-proton 
splitting of the electric polarizability is a iV" 1 effect (I), 
and we can calculate it consistently only to this order. 

In principle, one might try to perform a calculation 
which consistently takes into account the subleading 
pieces. The appropriate scheme would be the Kerman- 
Klein method [44j, but its application would involve a 
complicated fully quantum-mechanical calculation. 



VI. OTHER MODELS 

Techniques described in this paper are applicable to 
other model after straightforward modifications. In the 
Skyrme model, the described RPA method involves fluc- 
tuations of the meson fields which do not satisfy the non- 
linear constraint for the a and 7r field operators. This 
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linearization may be viewed as an approximation to the 
fully nonlinear dynamics. The RPA dynamics, obviously, 
involves mesons only, and the higher-derivative terms are 
manifest in the equations of motion for the fluctuations. 

In the case of the (partly bosonized) NJL model , 
the mesonic potential has the simple form i^ 2 (cr 2 + 7T 2 ). 
The sea quarks are present explicitly, and the number of 
quark equations is infinite. Standard methods of solv- 
ing these equations numerically may encounter problems 
for the case when the translational zero mode is excited, 
since extremely good accuracy is necessary in this case. 

VII. CONCLUSION 

We have presented the linear response method in 
hedgehog soliton models. We have shown that the 
method is consistent with the basic philosophy of these 
models, namely, the l/iV c -expansion, if its application 
is restricted to obtaining the leading- N c order of a given 
quantity. We have discussed many technical points which 
are encountered in practical calculations, especially the 
treatment of zero-modes, which create special problems. 
Appropriate equations of motion have been classified ac- 
cording to hedgehog symmetries, and derived explicitly 
for the model of Rcf . jD . Our method, after straightfor- 
ward modifications, is directly applicable to other hedge- 
hog models. A physical application of the approach is 
described in the preceding paper, (I), where we study 
the electromagnetic polarizabilities of the nucleon. 
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APPENDIX A: EQUATIONS OF MOTION FOR 
SMALL FLUCTUATIONS IN THE GRAND-SPIN 
BASIS 



We compose the basis of Dirac spinors with good K 
quantum numbers using the coupling scheme in which 
the isospin, I = |, and spin S = ^, are first coupled 
to a quantum number A, and then orbital angular mo- 
mentum, L, and A are coupled to K. Since there is no 
confusion concerning the value of K or K 3l we use the 
notation 

\L, A >= \K, (L, A(J =|,5= \)),K 3 > . (Al) 

States with parity V = (—) K (or V = — (— ) K ) are called 
normal (abnormal) parity states. The basis of Dirac 
spinors is 



Spinors X and Y are expressed in states ( A2 ) . The quark 



sources are decomposed into (L, A) components: 

3 L ' A =( . j <tl% ,)\L,A>, (A3) 

V 1(7 ■ r J F M / 

Tables Q - II] list the matrix elements which arise in deriv- 
ing the qua rk p arts of perturbation equations. It is clear 



from Table III that unless K — 0, the kinetic term mixes 
the A = and A = 1 components of the L = K states 
(normal parity case). Diagonalization is made through 
the substitution 



2K+1 V 2K + 1 



and similarly for the F-components, and the sources. 

The basis for the meson fluctuations is composed by 
coupling isospin to L. For a given value of K, the a and 
7r fluctuations can be expressed through functions 

a L (r)\K, (L, 0),K 3 >^ L {r)\K, (L, 1), K 3 >, (A5) 

Obviously, L = K for a, and L = K — 1,K,K+1 for tt, 
such that for a given K and K 3 

5<y = <j K {r)\K 1 (L,0),K 3 >, 

5tt= 7r L (r)\K,(L,l),K 3 > . (A6) 

L=K-1,K,K+1 

Using standard Racah algebra, it is straightforward to 



derive the general equations (3.4) for a given K pertur- 



bation. In the notation of this appendix, G^ XY y, etc., 



correspond to the X and Y spinors from Eq. (3.4) , 
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and G^x+y} = G °x + e ^ c -- The functions describ- 
ing meson fluctuation, a L , ir L , have the meaning of the 



Z-functions of Eq. (3.4) 



For normal parity equation we get 



K 

d r G a {X Y] = —G\ X y} + (g<?h - £ T ^)F{x,y) 



y v 2K+1 {X ' Y} 2K + 1 {X ' Y} ' 

K+l _ _ K K+l\ -a 



+ g{ \l WTl Fh a +Gh ^ )-3f,{x,y } , 
K + 1 

d r G {x ,Y} = G {X,Y} + (9<?h ~ E =F ^) F {X.Y} 



+ s M- 2K + 1 G {X , Y} + ^kTi ^ 

+ 9 ^ 2K + i Fh<jK ~ Gh^ K ^ X ) - 3f,{x,y}' 
drF a {x .y } = -^^F a {x , Y} + (9*h + £ ± u)G a {x , Y} 



+ ^2KTl {X ' Y} 2K + 1 

+ 9[ {Wk Gh(jK FhnK+1) + 

9rF b {x . Y} = ^- F ky} + fah + £ ± u)G\ xx] 



y V 2K + 1 {X > Y} 2K + 1 {XY}) 

+ 9i \l 2KTl Gh<jK + F ^ K ~^+iG,{X,Y}, 



(A7) 



+ 2A \2KTi h + V 2KTi <7h7rh ' T 2K + 1 * k * ) 

- gN c (F h G\ x+Y} + G h F h {x+Y} ) + £ ~\ 

(ti 1 _i_ 2 ,9 K{K+1) \ K 2 2 2 , , ,2x if 

+2A ^v^n^^ +<T/i " "V^n™" J 
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J VK(K + l) 2 K -! 

— — Tt h TT 



2K + 1 n V 2if + 1 

gN c (-F h G a {x+Y} - G h ff x+Y} ) + £+\ 



if + 1 K 



K+l r 
2K + Y 



The abnormal parity equations have the form 

if 



a r K-l,X 
Ur(jr {X,Y} 



l.oK-1,1 , / _,_ \£,K-1,1 



{X,Y} 

I „ K _ t>1 2y/K(K+l)„ K+h u 



G, 



2K + 1 W> 



2K + 1 



G 



{X,Y}> 



K + l 



2K+1 



Ghir K 



■ K-1,1 
■?F,{X,Y}> 



p. r K+l,l 
Ur ^{X,Y} 



if 



-G 



K+1,1 



{X,Y} 

2^/K(K + 
2K+1 



K+l.l 
{X,Y} 



1)^-1,1 _ 1 

u {x,y} 2K+l"^ x ^ 



-G 



K+1,1 



K 



G h 7T K - J 



2K+1 



■K+1,1 



d r F, 



K-1,1 
{X.Y} 



K + l 



-F, 



K-1,1 
{X,Y} 



F,{X,Y}' 



+ (ga h + £± to)G 



K-1,1 
{X,Y} 



1__ f k-i,i 2^K(K+l) ^ K+hl 



2IT + 1 W 

_K 



2K + 1 



F, 



{X,Y} 



) 



K + 1 F k , • 

2if + l l?r + Jg >W>' 



(9 F 



K+1,1 



^ pK+1,1 



, N^tK + 1,1 

±UJ ) U {X,Y} 



2y/K(K+l) 1 ^+1,1, 
ff7T/i( — — r '' 



2if 



1 ' pn-ia , 7-tJ\-h,i\ 
1 *{X,Y} 



A' 



•K+l.l 



2if + 1 



(pp- _l_ ^« if + K \2/ 2 , 2 ,2 , ,2\K 

[d r + -O r — ) 7T = A (<J h + TT h - V - UJ )TT 

- gN c [G h 
+ F h 



if + 1 _ f k-i,i 



2K +1 MW} T V 2K+1 <*+ y > 



if 



if + 



2if + 1^{x+y} -<- Y 2if + 1 



if ^K+1,1 




For the u> = 0, odd-7?. case, meson fluctuations van- 
ish, and appropriate equations have the form of Eqs. 
(A7,A9), with the meson fluctuations set to zero. The 
X and Y eq uatio ns can be combined to a single equation 
of the form (|3). 

In the case of an even- 1Z source which does not excite 



a zero mode (case Qo = in Sec. Ill B 1 ), we can set 
u> = 0, in the above equations. We ca n co mbine the X 
and Y equations, and obtain the form (3.8). If the zero 



mode is excit ed (Qo ^ 0) , we have to solve full equations 
(A.7- Aq), or (A9-AlCj), depending on parity. 
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For the special case of K = 0, G a = G 0,0 , G b = 0, etc., 



and only equations for the a components in Eqs. (A7) 
remain. Fields with negative (i.e. K — 1) superscripts, 
and equations for these fields are eliminated 



The boundary conditions in Eqs. (A7-A10) are such 
that the solutions are everywhere finite. At the origin, 
radial derivatives of S-wavc fields vanish, and the values 
of higher-L fields vanish. At r — > oo, the appropriate 
boundary conditions follow from solutions of the equa- 
tions in the asymptotic region. 



APPENDIX B: COLLECTIVE MATRIX 
ELEMENTS 

Suppose a space rotation, R, is described by Euler an- 
gles a, (3 and 7: 

R = e- iaJ *e- i/3J ye- i ' yJ \ (Bl) 



Then, matrix B from Sec. [II F is given by 



B = e n-r3/2 e i/3r 2 /2 e iar 3 /2_ (g 2 ) 

The matrix transforming from the body-fixed to lab 
frame, c a b, is defined as 



Cab 



\Tr[ Ta BT b B^] = Dl(a,P, 7 ), (B3) 



where the first (second) subscript in the Wigner D-matrix 
is connected to the spin (isospin) space. It follows that 
J2b Cab ^ b = ~ ^o- T ne s P m operator and the matrix c 
commute, [c a b, Jk] = 0. We also introduce a vector c 
defined as 

c=\Tr[T 3 BTB^]- Cfl =Dl . (B4) 

The collective baryon states with spin J, isospin I = J, 
and projections m and ^3 are 



|J = / ;TO) J 3 ) = y^±i£K (B5) 

In formulas below we do not display m or ^3 in labels of 
the states, and use notation \N) = |i;m, J3), and |A) = 
||;m, I3). The following useful formulas can be easily 
derived (no implicit summation over repeated indices): 

M*c M = i + (f-/i 2 )i} 2 , 5>/*)% = 1 ' (B6) 



from which follows that 

(N\{ Cfl )* Cf t\N) = i; any M , 

2 2 

<A|( C/1 )%|A> = i + 



3 ~ M J +1 ; |m| = |Z 3 | 



5 I -l ; \m\ + \h\ 

1 ; m = I 3 
1 ; m = -I3 



{N\(MA) = ^t^l!+]- : rn = h . (B7) 
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One also finds 

(N\((J^y^ + (^yj^\N) = -p 3 ; any M . (B8) 
One also derives 

(AT|co|A) = ^. (B9) 

For our analysis of the A states in hadronic loops in (I) 
the following formulas are important: 

(N\(c^r\N;m\I^(N;m\I^\N) = §, 
(N\(c ti y\A;m',I' 3 )(A;m>,l! s \c tl \N) = l (BIO) 

The sum of the a bov e formulas gives unity, in accordance 
to the sum rule (p36|). 



APPENDIX C: EQUALITY OF INERTIAL AND 
SOLITON MASSES 

For the case of translations, the inertia parameter, A4, 
is equal to the soliton mass, M so i. This result, required 
by Lorentz invariance, can be verified explicitly as fol- 
lows: Consider a boost in the z-direction, with small ve- 
locity v. The fields transform as 

<j) -> (f>h{r - vt), 

e -iet e - i£{t - Vz)+ ^ va 'q h (r - vt), (CI) 

which lead to the following shifts linear in the velocity: 

5<f> = -vtd z (f> h , 

5q h =ve- i ' t {iez+\a z -td z )q h . (C2) 

Using identities [h,z] — —ia z and {h,a z } = —2id z , we 
easily derive the equation 

[h - e](-ez + \ia z )qn = d z q h . (C3) 

After integrating by parts we get the expression for the 
energy shift of a moving soliton: 

58 = \v 2 {\T q + \T+ + N c e), (C4) 

where T q and are kinetic energies in the soliton, car- 
ried by the quarks and mesons, respectively. 

Next, we use a virial relation. Consider scale change of 
the radial coordinate, r — > sr. The soliton energy scales 
as £{s) —T q /s + Vq<f, + sT,p + s 3 V0, where V q< j, and V# are 
the quark-meson, and meson-meson interaction energies. 
Stationarity of the solution imposes d s £ \ s —\ = 0, which, 
together with the relation N c e — T q + V q< f,, leads to the 
virial relation 

M sol = \T q + \T+ + N c e. (C5) 



Comparing Eq. ( |C4J) and Eq. ( |C5| ) completes the proof 
that 5£ = \v 2 M so i. Using similar methods, one can 
show the equality of inertial and soliton masses in other 
models Em, also in non-local theories, such as the NJL 
model |471| . 
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APPENDIX D: BOSONIZATION AND PAULI 
BLOCKING OF THE DIRAC SEA 



In this section we return to the question whether the 
Dirac sea should be "Pauli-blocked" in our model. Effec- 
tive chiral models are believed to result from bosonizing 
QCD, which, of course, can only be done approximately. 
For definiteness, we discuss the issue of Pauli blocking in 
the framework of the partly-bosonized |4q] NJL model, 
but the result is more general. In presence of an external 
source, J, the action of the model is: 

Snjl — —iTr log[ifl — gU — J] — vac, 

gU = g(<r + ij 5 T ■ 7r), (Dl) 

A cut-off is understood, Tr denotes functional trace, 
and vac means the vacuum subtraction. For simplic- 
ity, we assume the nonlinear constraint cr 2 + 7T 2 = Fjr 2 . 
The source J may represent interactions with external 
probes (e.ff . electromagnetic) or result from cranking 



(Sec. HIE ). For definiteness, let us evaluate the mo- 
ment of inertia. In this case J = • r, and expanding 
the action to second order in A we obtain 

AS NJL = ±A 2 6 / dt, (D2) 



where the moment of inertia, 0, is given by 

6 = >c / Ps P ^-t't 3 , (D3) 



4 J 2n us — h lu — h 

where Sp denote the trace over space, spin and isospin, 
and h is the Dirac hamiltonian. The pol e st ructure and 



the contour of the u> integration in Eq. (D3) is given in 
Fig. (0). Note that the contour goes above the occupied 
valence state, as well as above all the negative-energy sea 



states. Performing the integration over us in Eq. (D3), 
we obtain the usual spectral expression for: 

e = i Nc £ mi^L, (D4 ) 

iGocc. 
jEunocc. 

where occ. denotes all occupied states, i.e. the valence as 
well as the sea states, and unocc. denotes the unoccupied 
positive energy states (see Fig. |l| for the meaning of 
labels). The expression under the sum is antisymmetric 
with respec t to exchanging i and j, therefore the sum 



as in Eq. (D4) over i and j belonging to the same set 
of indices vanishes. Using this trick we can replace the 
ranges of summation indices as follows: 

E = £'= £'+ £'= £'+ E 

?£occ. ■itEocc. iEval. iGsca i£val. z£sea 

jG unocc. jEaU j'Eall jEall jEall j Epos. en. 

(D5) 
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where the prime means the exclusion of i — j term, any 
denotes all states, and pos.en. denotes the positive energy 
states. According to Eq. (D5), the moment of inertia can 
be decomposed into the valence and sea parts: 



e = e 



val. 



e 



val 



l(*|r 3 |i)r 



5 ® sea — o^Yc 



^—^ P ■ — P 



iGsca 
jGpos.cn. 



(D6) 



Note that the "full" expression (D4 ) ob eys the Pauli ex- 
clusion principle, hence using Eq. (D5) we have broken 
the the original expression into two parts, each of which 
violates the Pauli principle. In fact, an analogous de- 
composition is used in the treatment of the relativistic 
fermion propagator in fermion matter |48f . Below we ex- 
plain why this is useful. Firstly, the expression for Q va i 
corresponds to our quark part of the moment of inertia 
calculated in Sec. Ill Ej . Secondly, the sea part of the 
moment of inertia ca n b e simply approximated only if it 
is written as in Eq. (D6). Indeed, we can write down 



e., 



~ A N C / — Sp 

4 2ir uj — h u) 



1 



-r 3 - 



(D7) 



where the contour of integration is given in Fig. H This 
contour can be Wick-rotated without picking up any pole 
contributions, and we obtain an expression in Euclidean 
space. We can then perform standard gradient expansion 
methods p9| , pO|j5l|Jl7| ] to rewrite 8 sea as an integral over 
the classical pion field. The first term, with no deriva- 
tives, is just our expr essio n for the pion part of the mo- 
ment of inertia, Eq. ( |3.2l| ). Furthermore, this term does 
not depend on the NJL cut-off, since the normalization 
factor is the same as in the pion wave function normal- 
ization p5| . Further terms in the gradient expansion do 
depend on the cut-off. If we tried t o p erform the Wick 
rotation on the original expression (D3), we would pick 



up a pole contribution fr om the occupied valence level, 
and our final expression (D6) would also follow. 
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FIG. 1. Contour of integration, C, for the total (sea- and valence-quark) contribution: C cannot be Wick-rotated without 
picking up the valence quark contribution. Notation for various labels used in the text is visualized. 

FIG. 2. Contour of integration for the sea-quark contribution: C can be Wick-rotated to the contour C. Upon bosonization, 
the sea-quark effects can be described by mesonic degrees of freedom. 
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TABLE I. Matrix elements of r ■ r 
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TABLE II. Matrix elements of a ■ r 
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TABLE III. Matrix elements of a ■ L 
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